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Let Ibea space of functions, say X C C(K), K locally compact Hausdorff, let leX* be an integral 
on X and let M* C X* be a given subspace of "simple" functionals, then it is desired to obtain an /eM* 
for given rc, /eM* C M*; M* being a suitable n dimensional subspace determined so that!/ — / 
annihilates a given finite dimensional subspace X x C X. In this general context, the abstract analysis 
of numerical integration is developed and certain specific applications are made. 
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1. Introduction 

In the development of numerical integration formu- 
las, one is given a function space X, say X C C{K), 
where K is some locally compact Hausdorff space, 
and an integral leX*. It is then desired to find an 
leM* C A'*, where M* is a subspace generated by a 
family of "simple" integrals {x*} such that / can be 
said to approximate / in some given sense. 

In ordinary quadrature, the family {x*} are given 
by x*f=fix a )Yj'eC(K) where x a eK. An element 
leM* is said to approximate the given integral / if 
1 — leXi where X^ is the subspace in X* of all anni- 
hillators of a given subspace X^ C X. 

The following paper summarizes the essential mathe- 
matical structure of numerical integration methods 
with the purpose of extending systematically the class 
of such methods. 

2. General Analysis 

Let X C C(K) and let X () C X be a subspace 3 dim 
X () = n and X () is generated by a given set of n "basis" 
functions {fk}- In X* take M * a subspace generated 
by a given family {x*} of linearly independent func- 
tionals with the properties: M* is total over X and 
{x*} H Xi is void where X& = {1; leX*, Z/=0,/tf }. 
If M* C M* designates a subspace generated by n of 
the x*, say {x a i}i = l ... /i so that dim M* = n, then 
one has X* =X^ -f-Af* since the deficiency of X^ = dim 
X =n. Also M* is total over X as is readily seen. 

If leX* is a given integral, then } T eM* bI = T+ l\ 
where l^eX^ and } such an I for each subspace M* 
as defined above. For any such given subspace M* 
} li eM* i = l ... reBlif k =8 ik i 9 k=l...n d ik =l 
i = k, 8tk=0 &4= k. This follows from M* being total 
over X(). 



Thus one can write: / = V /(//)//+ /^, l\eX^ and 

i 

li = ^ 7ijX* where x* = x* • • 
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(1) 



The functional / therefore depends on the selected set 
{x*(\ , that is on M* n on the If and on the coefficients 
y, } i,j = l . . . n. 

The relations V y^x *fk = 8iki, k=] ...n can be 

j 

put in matrix form: 110=/ where / is the identity 

matrix, r=(y, y ) r {) =(x*f\). Now V {) is nonsingular 

for any given M* since M* total over X {) implies that 

^ TkX*fk = j= 1 ... /i is possible only when ta = 
A-= i 

/c = 1 . . . n. Hence the matrix r= I^ 1 is determined 
once To is known. 

If the {fk} are such that the //a are readily evaluated 
then the problem of determining / is that of specifying 
M* or a corresponding basis {#*/}. Under the con- 
ditions imposed thus far with dim X () = dim M* = n 
/ will not in general be unique without additional 
conditions on X and/or the M* 



3. Specific Cases 

In certain quadrature procedures, such as the New- 
ton-Cotes formulas, M* is specified a priori by 
choosing the x* } . The subspace Xq in these 
cases is such that the If are easily calculated and thus / 
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is determined. For instance, the Parabolic Rule has 
I— dx, where X is the space over the Reals gen- 

JXi 

erated by thef k = x k ~\ k = 1,2,3 and x*/=J{xi) for given 
real numbers X\, X2= 2 (X1 + X3), x 3 . 

In Gaussian quadrature methods, the following 
additional conditions limit the choice of M*. Consider 
a subspace X\ C X, where X\ D X and X\ has genera- 
tors {fk}, k = 1 . . . p, p > n, fk =fk for /c = 1 . . . n. 
Now if T is a linear transformation 011I3 f k = Tf k _ l 
= T k ~ l f[. Then T induces a linear transformation 
Q on J* given by Qlf=lf=lTf, V/eY*, /e¥. 
Consequently x*f k = (Qx*)f k -i = {Q k ~ x x^)f[ and one can 
write: 



(2) 



to = S (2 *M < ? fc " 1 ^/' for A= 1 ... p. 



Since the matrix r = (x*f k ) = (Q k ~ 1 x*fi) for 7, 
A: =1 ... 71 is nonsingular, the /i-tuple (Q n x*fi) 
j = 1 . . . n corresponding to 



(/» + >=2(X(/b'oj<2 n *;/. 

j \ i 
is linearly dependent on the rows of To, i.e. 

0"*//i=2 P^xffi j=l ■ ■ ■ n. (3) 

/V=l 



This gives the relationship 



(4) 



C" - 2 W-' )%*/;' = 7=1 

^ 1 = 1 

will imply that 



Now, if for Mf lL = {f; feX u s*/=0, ;=1 . . . n) 
one also assumes: TM^CM*^ that is, xff=0 implies 
x?Tf=0 = Qx*f for 7=1 . . . n then this condition, 
plus the statement: 



(5) 



<? n+ «-2)3*<2* + «-i *f/;=o y=i . . . »,V9 

(6) 
Taking g = 0, 1 . . . ft— 1 one gets the ft equations 

In Gaussian quadrature, the stronger condition 
/ — /^A^-CA^- is imposed. Thus, if l — T*X\, if the 
,/JJ+q are such that the lf k+q are readily obtained and if 
the matrix L = (lf k+q ), k=l . . . n, q = . . . ft— 1 
is nonsingular then the /3a- can be obtained uniquely. 



Given the fik, solution for the xff[ 7= 1 ... ft follows 
from the relations (5) when they are solvable. Deter- 
mination of the functionals x* 7 = 1 . . . n proceeds 
from the values of the xff\ when the xf are sufficiently 
simple, for example, a one parameter family like those 
in the examples below. 

Now, L = (lf k+q ) is nonsingular if and only if ^ 
scalars 



e, + 0,A = l . . . nBl[^€ k f k+q =0 

q=0 . . . Ti-1. 
Since these last relations can be written as: 



l (2 e :Jk I*-i+f/ij = Z (T« (2 e tP-'J/^O (8) 

one has, in this general context, an orthogonality con- 
dition stating that no ft — 1 degree "polynomial" in T 
is orthogonal in this sense to all ft — 1 degree poly- 
nomials in T. The ft-equations (7) can also be written 
as: 

/(Wr»-]T ^-^/iW q = . . . ft-1 

(9) 

which states that the n degree polynomial T n — 
V ftkT k ~ l is orthogonal, in this sense, to all ft — 1 

A = l 

degree polynomials in T. 

Thus if Xo, Xi, M% satisfy the conditions imposed, 
it w_as seen in the general section above thatJ/eM* 9 
1 — IgXq. The preceding discussion shows that if 
l-UXi then (9) follows. Conversely, if 3- p k 9 (9) 
holds one shows directly l — leX{ using (5) and (6). 
Consequently one has the following statement ana- 
logous to the ordinary Gaussian quadrature theorem 
[1] which is usually stated in terms of classical orthog- 
onal polynomials: If X , X u Mf x satisfy all the condi- 
tions imposed above then} TeM* 9 / — leX{ if and 

n 

only if J- an n degree "polynomial" T n -^ p k T k ~ l 

k=l 

orthogonal as in (9) to all polynomials in T of degree 
<ft. 

4. Applications 

In ordinary Gaussian quadrature, the above spe- 
cializes to the case where X ={x i ~ 1 ; i=l . . . ft}, 
Xi — {x i ~ 1 ;i=l . . . 2ft}, and the {x*} are defined by 
x%f=f(xa). Then the other Gaussian conditions are 
seen to hold when T is 9(Tf)(x)=f{x)x. 



An interesting nonclassical application is to /= I dx 

on the space X=C (K)={f; feC(K), K=[0, l]/(0)=/(l) 
= 0}. The natural basis functions here are sin rrkx 
or the equivalent functions f k (x) = sin irx (cos 7rx) k ~ l . 
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Take Xq, X\, with this basis and x* defined by 
Xaf=f{oc a )- > x a e[0, 1]. Take T?Tf(x) = f(x) cos ttx 
then TMt±CMti and one gets: 

7=0.39 (xf + x$) for n = 2 

1= 0.205 (xf + x$) + 0.247 (x 2 * + xf) for n = 4 

7= 0. 18502 (*? + x$) + 0. 12306 (x 3 * + * 4 *) for n = 6. 

+ 0.16882 (* 5 * + *e) 

Another application is to 1=\ k(x)dx for some 

Jo 
weight function k(x) e.g., k{x)= e~ x2 . Here K= [0, <»] 

Z= ]/; feC(K), //exists, /fr)«ft exists 

and is easily obtained for arbitrary *eA f• 
Take/,;(x)=A;* /c - 1 , let xff=k{&) \ 3 f(x)dx. For x 
= Xj-i<ijj<Xj and take T given by Tf(x) = xf(x) 

H- r ./UW/ then Jl^CM^ and one gets: 
Jo 

-0.922 | fib + 0.195 J * dx 

for n = 2. 



where s = 0, jt, = 0.754, x 2 = 1.734 



fj-i r« f^ 

/ = 0.9798 d* + 0.6188 dx + 0.1386 <fo 

J^O Jj"l Jj'2 



Jx 3 



+ 0.0051 dx 

Jx 3 

for n = 4> 



where * = 0, xi = 0.4239 * 2 = 1 .014 

^3 = 1.792 * 4 = 2.640 

f xi r xz r xz 

= 0.99200 </*• + 0.81425 r/*-h0.41032 dx 

J.r„ Jx, Jx 2 

+ 0.09445 \ dx 

C XS f xe 

+ 0.00702 dx + 0.00008 <2x 

J J"4 J .*'.-> 



for n =6 

where i o =0i,- 0.27779 x 2 = 0.68991 

^3=1.2090 x 4 = 1.8138 x 5 = 2.5104 
x 6 = 3.3556. 

Greater detail on the development of such formulas 
and their utility will be given in a subsequent paper. 
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